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SUMMARY 


A method for ccanputlng the three inertial components of a previously 
unknown lunar landing site is presented. The conqputation of the landing 
site vector csm be sin5>ly computed onboard the CSM while it is in orbit 
about the moon prior to separation for descent. The computation requires 
that two unit vectors, separated by a few minutes in time, be measured with 
an optical instrument to the landing site and with respect to the IMQ. The 
method is also extended to include the computation of the expected errors 
in the landing site i>osition. These errors are dependent on the expected 
CSM errors in position and velocity as well as on the error in the sighting 
instrument. An example is presented to illustrate the computation. 


INTRODUCTION 


For the Apollo lunar landing mission the landing site coordinates 
are expected to be selected well in advance in ord^ to insure adherence 
to the many operational constraints imposed on the mission. However, the 
operational flexibility exisns in the onboard guidance computer program 
to allow the navigator, while in liuiar orbit, to select a new landing site 
or to update the preselected site. Due to the large uncertainties that 
exist in the current estimates of the lunar radius, it is anticipated that 
an iq>dece of the landing site altitude alone would substantially increase 
mission success. 

This note will present a method for determining the landing site vector 
in inertial coordinates, during the orbit phase of the Apollo I"nar Landing 
Mission, using the onboard optics as they are presently planned. 

This note will also include as part of the method the derivation of a 
set of equations from which the e:q)ected errors of the landing site vector 
may be calculated. These errors may be conputed either on the earth, or 
onboard in the event that the expected errors are found to be a necessary 
input fcxr the LQ4 descent guidance scheme. An example case is presented 
for illustration. 


Determination of the Landing Site Vector 

At a time, t^, prior to the lunar landing, and while the CSM and LEM 

are both in orbit about the moon, the navigator sights the optical instru- 
ment on a distinguishable future which is located on the lunar surface in 
the general area where ^e landing is to occur. When the mark button is 
pressed, a unit vector, ^ , from the spacecraft to the laru^ing site is 
obtsdned. This unit ve^or is determined to within the sighting accuracy 
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by the two angles, C6^and , which are measured by the sighting 
instrument with respect to the Inertial Measurement Unit (IMU); see figure 

A 

1. The unit vector, ^ , is given by, 

^ ^ ^ t/>S ( 1 ) 

A ^ \ 

The unit vectors ( , j , K , are along the selenocentric (inertial) 
coordinate sixes. As presented here, it is assumed that the IMU is sdlgned 
along ^ f J f f so that oC and ^ refer to right siscension and declination. 
There is no loss of genersJ.ity in naking this assumption. 

The navigator continues to track the landing site through a period 
of a few minutes. At time t^, the navigator again presses the mark button, 
determining in the same msmner the unit vector, 

- < CH>sS, COSaC.+j tOS St ^ 

Since the moon is rotating very slowly, the following relation holds 
to a good approximation, 

4 - !Cc-* ^ f.f. (3) 


as can be seen from figure (l). 


In order to determine the landing site vector, St it is necessary 
that the magnitude of the vector ^ or ei be known. There is no way of 
deter mining the magnitudes p or P| by measurement. However, knowing 
Pjj and (>, , and the position vectors JT^ and or may be 

calculated. This may be shown as follows: Fr(mi equation (3), 




(U) 
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Take the dot product of (U) with ^ , 


Take the dot product of ( 4 ) with ^ , 

A A 

multiply (5) t>y ^ jo' p result to (6), 

or by simplification, 

^ ^ C,. C<as}«5; 

{j> - — r-^— 

SA4J 

/ A A 

where, dos p = 

- / - Cfi-’Pf 


With ^ determined by ( 7 )# the landing site vector xoBy be determined by 
combining (7) and (3) 




S/aJ 
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In summary > that is required to compute the landing site vector 

are the spacecraft position vectors at t^ auid t^ amd the measured 
quantities , 5^ ,06, , and * 


Determination of the Estimation 
Errors in the Landing Site Vector 


Let the actual position of the lamding site he given by. 


A 





( 9 ) 


where and ^0 are the actual position vectors of the spacecraft and 
the lamding site with respect to the spacecraft at time t^. Now let 
equation (3) represent the similar relation between the estimated landing 
site position, A , the estimated position vector, , and the estimated 
vector, ^0 , from spacecraft to landing site. 


That is. 


4 = fo 


where 



(10) 


The error in the estimate of ^ is n<nr found by subtracting (9) from (lO)y 

A- rf 


$4 - -f- ^£0 


( 11 ) 
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The covariance matrix of estimation errors in the landing site vector 
may he found by taking the expected value of (ll) times its transpose. 




* t (sr^ 


( 12 ) 


The term !Zjis the covariance matrix, , of estimation errors 


in spacecraft position at t^. Define, 




(13) 




(ii^) 


Then (l2) becomes. 




Co^o 


T 

( 1^0 


It is shown in Appendix A that equation (15) reduces to. 


(15) 


-h A//^A/^ 


(16) 
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where is the covariance matrix of estimation errors in spacecraft 
position at time t^. The other quantities in (l6) are: 



where CT is in the standard deviation of the sighting instrument. 








( 20 ) 


where is the cross correlation between velocity and position 


estimation errors at time t^. The matrices and ^ 2 . sulanatrices 

of the six by six transition matrix for Keplerian luotion; for exanq>le^ 
see reference 1, 


The matrix may be conputed from 




where ^^is the covariance matrix of estimation errors in spacecraft 
velocity at time t^. 


( 21 ) 
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In summary, all that is required to coinpute the estimation errors 
in the landing site vector is the covariance matrix of estimation errors 
at t^; the measured quantities oCa standard 

deviation (T of the sighting instrument. 

A Numerical Example 


The equations to compute the landing site vector and its errors were 
programmed on a digital computer. Inputs to the program were taken from 
the results of a completely independent program. The following quantities 
were given to the program: 

Ephemeris Time (ET) = I969 yr, 260 day, 77 hrs, 10 mins, 12.2 secs 
t = 0.703125 hrs from ET 

V(? 

t^ = 0.765625 hrs frcm ET 

At = (-93^.952, 370.206, 183.861) n.mi. 

(3^0 = 21.3^39° 

^ = 14.4697° 

At t^, X, = (-837.079, 527.752, 252.152) n.mi. 

= -85.3062° 

= -4 o. 13 i 4° 

With these quantities given, the inertial components of the landmark were 
found from equation (8) to he. 


= (-5.02893, 2.50418, .936535) X 10^ ft. 
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In order to compute the estimation errors in the landing site vector, 
additional statistical quantities were also input as follows: 


Er 


o 



s 

96819 

.840991 

-. 0841973 I 



4. 55638 

-. 4947417 1 


(SYM) 


5.475946 J 


^'.751570 

.0559336 

.2127881 



1.491713 

-. 32559 ^ 


(sym) 


8. 173442 J 


5 2 
X lO^ft 



2 

sec 


The root-mean- square position and velocity errors computed from the traces 

of Er and Ev^ are found to be, 
o o ' 


RMSPOS(t^) = 1095 ft 


RMSVEL(t^) = 3.23 ft/sec 


.580603 

. 3 ^^ 95^3 

. 7^310 


2.06510 -Al79^1 p ^2 
3.8873^ -6.429013 / X 10 — 
-6.83694 17.69110 J 


The sighting accuracy of the instrument was chosen to be 


^ = .001 rad 


Using these inputs, the landing site covariance matrix was found from 
equation ( 16 ) to be. 




4.59429 2.88405 

5.75483 


-1.90845 e p 
- 1.48731 X lO^ft 
6.28671 


The root-mean-square landing site error is confuted from the trace of E_ 

L 

and is found to be. 


RMSLS = 1290.0 ft 
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If E is transformed into the solenographic polar coordinate system as in 
L 

reference 2, the root-mean-square errors in latitude, longitude, and vertical 
directions are given by. 


RMSIAT = 85^^. 6 ft 
BMSLON = 800.2 ft 
RMSVERT = ft 

In reference 3 it is seen that the predicted uncertainties of known lunar 
landmarks vary from 220 meters to 910 meters in the horizontal and from 
730 meters to 800 meters in the vertical. The uncertainties in the 
vertical are possibly even larger. Based on results from the Ranger series, 
the Jet Propulsion Lab has recently decreased its estimate of the lunar 
radius from 1738 »0 KM to 1735 KM. 

The results of this example show that these uncertainties can be 
significantly Improved during the Apollo Lunar Landing Mission. 

CONCLUDING REMARKS 

The analytic description of a method for determining an unknown lunar 
landing site vector during the lunar orbit phase of the Apollo mission 
is presented. The method requires the measurement using an onboard optical 
instrument of two xmit vectors to the landing site with respec'' to the IMU. 
The method includes the computation of the landing site errors, which are 
dependent upon the CSM position and velocity errors, as well as the error 
in the sighting instrument. An exanrpie was presented to illustrate the 
computation. The results of this example show that by the use of onboard 
sightings on the landing site, the errors in the components of the landing 
site vector are determined to be substantially smaller than predicted 
luicertainties for lunar landmarks. 
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API^IX A 


The Covariance Matrix of Estimation 
Errors in the Landing Site Vector 


CcHisidcr the equation (15) for the covariance matrix of estimation 
errors in the landing site vector. 




• 4 - S” 

(a 




(Al) 


where > and y o ^ are given ty (13) €Uid (ik). In order to evaluate 
(<» ^ 

and it is expedient to find an esqpression fcr that is 

easily differentiated. Using (7), 


L J 


A 


A. A 


A A . 2 


using (U>s4 ^ 


e. = e. [L- (e.- e.)e'f [i -(r.- e.f/"c, 




or hy defining. 




A ^ T r •T' ^ ^ T 'I 

6^ '' fi J 


A A I 2 


^ ^ ^ J 


(A2) 





(A-2) 


The «q;ire88ion for bec(»ie8^ 




f - M £. 


(A3) 

Hbv 

the fiT8t order deviatlonB to hecone. 




-f- A/Xe 

(Alt) 

where 

5 c, =. Sr, -Sr^ 


But 

ll 

S^i ^ 

(A5) 

80^ 


(a6 ) 

Hbv define the vectore. 




- €^a 

^ ^ c. 

<^ 0 ^o ~ 



-=^ <L 

s?s, - 

, -£4 = •^ s 

tpcC, 


1 

1 

may be written. 




^ (t> 

+ MSSl 

(AT) 

where 

A/ = r-??, ,-s2 

(A8) 



(A-3) 


and. 



(A9) 


Nov consider equation (l4)^ using (AT) and 


Sc = S£,-Sr^ , 


- A/ B(S^Sr; ) 

- M£(^r, SrJ) 


(AlO) 


The first term of (AlO) vanishes if it assumed that the measurement errors 
i?nd trajectory errors are uncorrelated. 


e(s^sr7)^o 

(All) 



as shown in reference 1, the second term in (AlO) heccmies. 


(A12) 


S(^r,St7) = je,£(Sr.SrJ)^ Scj) 




(A13) 



where Br r is the cross correlation between velocity and position estimation 
errors . 

Equation (AlO) finally becomes, 

(Alu; 

Now consider equation (13), using (A7) 

If, = 

- £[( A/Ki +M ^Sr)(S^ V 

+ a/e(s^SsJ)m^ 

(A15; 

The third and fourth terms vanish by invoking the assumption that measurement 
and instrment errors are uncorrelated, that is, 

£(_S^ %$£) -- 8 (Sr, sr)- £(s-r = o 


Define the matrix 

- E(s:^s^y 


(AlC) 


which can be shown using reference 2 to be 

f^-J^D//f<f(tjXc^S„, !j ^ec^S, ) 


(A17) 



(A-5) 


Using 5-c - Sr,- fr, eq.\iation (A15) becomes, 

M[S C^r, fj) 

- SC^r^ $rrj~ ^r/jj M'" ^ A/^a/ 


All of the terms in the last equation have been previously defined, so 
finally equation (13) becomes 

Zf. -- - 5s - (^r: f--) 

M’’ -f zi/^vi/’" 

(A18) 


Using the definition (20) for the matrix D, along i«ith equations (A14) and 
(A18) the equation (A15) becomes 


=^f;- 




(A19) 




Figure la.- The geometry that defines the relations betwe^ landing 
site vector, ^ ; the measured unit vectors and ^ 
and the spacecraft position vectors at t and t^. 





